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Abstract 

^3 ■ 

We analyze approximation algorithms for several variants of the traveling sales- 
man problem with multiple objective functions. First, we consider the symmetric 
. TSP (STSP) with 7-triangle inequality. For this problem, we present a determin- 

' istic polynomial-time algorithm that achieves an approximation ratio of min{l + 

^0 ' 7, 2^:^^327+1 ^ ^ ^ randomized approximation algorithm that achieves a ratio 

' of ^^2^2^+1 ~^ ^' "'^^ particular, we obtain a. 2 + e approximation for multi-criteria 

c/3 ! metric STSP. 

I ^ I ' Then we show that multi-criteria cycle cover problems admit fully polynomial-time 

randomized approximation schemes. Based on these schemes, we present randomized 
CO ' approximation algorithms for STSP with 7-triangle inequality (ratio i^^^24j2 + e) , 

(3 I asymmetric TSP (ATSP) with 7-triangle inequality (ratio 5 -I- ^J^^i + e), STSP with 

' weights one and two (ratio 4/3) and ATSP with weights one and two (ratio 3/2). 

o ■ 
\o ■ 

^ ; 1 Introduction 

O ■ 

I In many practical optimization problems, there is not only one single objective function 

O ' to measure the quality of a solution, but there are several such functions. Consider for 

^ , instance buying a car: We (probably) want to buy a cheap car that is fast and has a good 

gas mileage. How do we decide which car is the best one for us? Of course, with respect to 



X 



^ • any single criterion, making the decision is easy. But with multiple criteria involved, there 



is no natural notion of a best choice. The aim of multi- criteria optimization (also called 
multi-objective optimization or Pareto optimization) is to cope with this problem. To 
transfer the concept of a best choice to multi-criteria optimization, the notion of Pareto 
curves was introduced (cf. Section 1.1 and Ehrgott [12]). A Pareto curve is a set of 
solutions that can be considered optimal. 

However, for most optimization problems, Pareto curves cannot be computed effi- 
ciently. Thus, we have to be content with approximations to them. 
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Onhne Algorithms (WAOA 2006) [20]. 
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The traveling salesman problem (TSP) is one of the best-known combinatorial opti- 
mization problems [16, 18]. An instance of the TSP is a complete graph with edge weights, 
and the aim is to find a Hamiltonian cycle (also called a tour) of minimum weight. Since 
the TSP is NP-hard [14], we cannot hope to always find an optimal tour efficiently. For 
practical purposes, however, it is often sufficient to obtain a tour that is close to optimal. 
In such cases, we require approximation algorithms, i.e., polynomial-time algorithms that 
compute such near-optimal tours. 

While the approximability of several variants of the single-criterion TSP has been 
studied extensively in the past decades, not much is known about the approximability 
of multi-criteria TSP. The classical TSP is about a traveling salesman who has to visit 
a certain number of cities and return back home in a shortest tour. "Real" saleswomen 
and salesmen do not face such a simple situation. Instead, while arranging their tours, 
they have to bear in mind several objectives that are to be optimized. For instance, the 
distance traveled and the travel time should be minimized while the journey should be as 
cheap as possible. This gives rise to multi-criteria TSP, for which we design approximation 
algorithms in this paper. 

1.1 Preliminciries 

Graphs and Optimization Problems. Let G = (V, E) be a graph (directed or undi- 
rected) with edge weights w : E ^ 1^. We define the weight of a subgraph G' = {V',E') 
of G or a subset E' of the edges of G as the sum of the weights of its edges: w{G') = 
w{E') = J2e&E' ^(e)- For keN,we define [k] = {1,2,..., k}. 

TSP in general is the following optimization problem: Given a graph with edge weights, 
find a Hamiltonian cycle, i. e., a cycle that visits every vertex of the graph exactly once, of 
minimum weight. In case of undirected graphs, we speak of the symmetric TSP (STSP), 
while in case of directed graphs, we refer to the problem as the asymmetric TSP (ATSP). 

An instance of A -STSP is an undirected complete graph G = {V,E) with edge 
weights w : E ^ N that fulfill triangle inequality, i.e., w{{u,v}) < w{{u,x}) + w{{x,v}) 
for all distinct vertices u,v,x G V. 

For 7 G [i, 1], A (7) -STSP is the restriction of A -STSP to instances that satisfy 7- 
strengthened triangle inequality, i. e., w{{u, v}) < 7- {w{{u, x}) +w{{x, v})) for all distinct 
vertices u,v,x. 

STSP(1, 2) is the special case of A -STSP where only one and two are allowed as edge 
weights, i. e., w : E ^ {1,2}. 

A -ATSP, A(7)-ATSP, and ATSP(1,2) are defined like their undirected coun- 
terparts A -STSP, A(7)-STSP, and STSP(1,2), respectively, except that the graphs are 
directed. 

Note that for 7 = 1, A(7)-STSP and A(7)-ATSP become A -STSP and A -ATSP, 
respectively. As 7 gets smaller, the edge weights become more and more structured. For 
7 = 1/2, all edge weights arc equal. The 7-strengthened triangle inequality can also be 
considered as a data-dependent bound [7]: Given an instance of metric TSP, we compute 
the minimum 7 such that the instance fulfills 7-strengthened triangle inequality. If 7 < 1, 
then we obtain a better performance guarantee for our approximate solution than with 
triangle inequality alone. 

A cycle cover of a graph G = (V, E) is a subgraph (V, G) that consists solely of cycles 
such that every vertex v e V is part of exactly one cycle. In most cases, we refer to a 
cycle cover as the set G of its edges. Hamiltonian cycles are cycle covers that consist of 
only a single cycle. 
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The problem of computing cycle covers of minimum weight in undirected graphs is 
called sec. The directed version of the problem is called ACC. 

Multi-Criteria Optimization. A A;-criteria optimization problem consists of a set / of 
instances, a set sol(a;) of feasible solutions for every instance x e I, k objective functions 

wi, . . . ,wig, each mapping pairs oi x £ I and y G sol(x) to N, and k types indicating 
whether Wi should be minimized or maximized. We refer to Ehrgott and Gandibleux [12, 
13] for surveys on multi-criteria optimization problems. Throughout this paper, we restrict 
ourselves to problems where all objective functions should be minimized. Furthermore, 
we assume that the number k of criteria is fixed. The running-times of our algorithms are 
exponential in k. But since k is typically a small number, this does not cause any harm. 

The optimization problems defined in Section 1.1 are generalized to their multi-criteria 
counterparts in the obvious way: We have k objective functions wi,... ,Wk, each induced 
by edge weight functions (to which we also refer as wi, . . . ,Wk) as described. If we have 
additional restrictions on the edge weights, like the triangle inequality, every edge weight 
function is assumed to fulfill them. 

In general, the different objective functions are in conflict with each other, i.e., it is 
impossible to minimize all of them simultaneously. Therefore, the notion of Pareto curves 
has been introduced. For the following definitions, let 11 be a A;-criteria optimization 
problem as defined above. 

A set V{x) C sol(a;) is called a Pareto curve of x if for all solutions z G sol(a;), there 
exists a solution y G V{x) with Wi{x,y) < Wi{x,z) for all i G [k]. 

A Pareto curve contains all solutions that might be considered optimal. If there are 
two solutions y and z with Wi{x,y) = Wi{x,z) for all i G [k], then it suffices to put one 
of them into V{x). For completeness, let us mention that Pareto curves are not unique 
in general: In our definition, it is not forbidden to include dominated solutions in V{x) 
(a solution y is dominated if there exists a z with Wi{x,z) < Wi{x,y) for all i G [k] and 
Wi{x, z) < Wi{x, y) for some i G [k], i. e., z is strictly better than y). 

For the majority of multi-criteria problems, computing Pareto curves is hard for two 
reasons: First, many two-criteria problems allow for a reduction from the knapsack prob- 
lem. Second, Pareto curves are often of exponential size. Therefore, we have to be content 
with approximate Pareto curves. Let /3 > 1, and let a:; G / and V^'^^{x) C sol(x). The set 
•papx^^-j -g called a /3-approximate Pareto curve for x if, for every z G sol(a:), there 
exists a y G V^^^{x) with Wi{x,y) < P ■ Wi{x, z) for all f G [k]. 

A 1-approximate Pareto curve is a Pareto curve. For completeness, let us mention that 
if n is a maximization problem (or an objective Wi for some i G [k] should be maximized), 
then the condition is Wi{x, z) < P ■ Wi{x, y). 

While Pareto curves itself are often of exponential size, it is known that (1 + e)- 
approximate Pareto curves of size polynomial in the input size and 1/e exist [23]. (The 
technical restriction is that the objective functions are restricted to assume values of at 
most 2Pi\^\) for X G / and some polynomial p.) 

The above definition leads immediately to the notion of an approximation algorithm 
for multi-criteria optimization problems: Let (3 > 1. A /3-approximation algorithm for 
n is an algorithm tliat, for every input x G /, computes a /3-approximate Pareto curve for 
X in time polynomial in the size |x| of x. 

A randomized /3-approximation algorithm for 11 is a polynomial-time algorithm 
that, for every input x £ I, computes a set 'P'^p^(x) C sol(x) such that V^'^^lx) is a 
/^-approximate Pareto curve for x with a probability of at least 1/2. 

By executing a randomized approximation algorithm i times, we obtain a /3-appro- 
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Variant 


Ratio 


Reference 


A-STSP 


3/2 


Christofides [10] 


A(7)-STSP 


mill -f „ ^'^T — T . „^~J7 \ 
1 37^-27+1 ' 3-37 J 


Bockenhauer et al. [8] 


STSP(1,2) 


8/7 


Bernian, Karpinski [5] 


A-ATSP 


0.842 • log n 


Kaplan et al. [17] 


A(7) -ATSP 


min{2z^,iZ^} 


Blaser et al. [7]; Chandran and Ram [9] 


ATSP(1,2) 


5/4 


Blaser [6] 



Table 1: Approximability of single-criterion TSP. 



ximate Pareto curve with a probability of at least 1 — 2~ , i.e., the failure probability 
tends exponentially to zero: We take the union of all sets of solutions computed in the £ 
iterations and throw away all solutions that are dominated by solutions in the union. 

Given the notion of (randomized) approximation algorithms, we can define approxi- 
mation schemes. A fully polynomial-time approximation scheme (FPTAS) for 11 
is an algorithm that, on input a; € / and e > 0, computes a (1 -|- e)-approximate Pareto 
curve in time polynomial in the size of x and 1/e. 

A fully polynomial-time randomized approximation scheme (FPRAS) for 11 
is a randomized approximation algorithm that, on input x G I and e > 0, computes a 
(1 + e)-approximate Pareto curve in time polynomial in the size of x and l/e. 

Finally, we define the notion of a randomized exact algorithm: A randomized exact 
algorithm for 11 is an algorithm that, on input x, computes a Pareto curve of x in time 
polynomial in the size of x with a probability of at least 1/2. 

An optimization problem 11 is said to be polynomially bounded if there exists a poly- 
nomial p such that the following holds for every objective function Wi of 11: For every 
instance x and every feasible solution y for x, Wi{x, y) < p{\x\) for all z G [k]. Analogously 
to the fact that a polynomially bounded single-criterion problem that admits an FPTAS 
can be solved exactly in polynomial time (cf. Ausicllo et al. [3, Theorem 3.15]), random- 
ized exact algorithms exist for polynomially bounded multi-criteria optimization problems 
that admit an FPRAS. 

1.2 Previous Results 

The approximability of single-criterion TSP has been studied intensively in the past. Ta- 
ble 1 shows the currently best approximation ratios of the variants for which the multi- 
criteria counterparts are considered in this paper. 

While single-criterion optimization problems and their approximation properties have 
been the subject of a considerable amount of research (cf. Ausiello et al. [3] for a survey), 
not much is known about the approximability of multi-criteria optimization problems. 

Papadimitriou and Yannakakis [23], by applying results of Barahona and Pulley- 
blank [4], Mulmuley et al. [21], and themselves [22], showed that there exist FPTASs for 
multi-criteria minimum-weight spanning trees and the multi-criteria shortest path prob- 
lem and an FPRAS (more precisely, a fully polynomial RNC scheme) for the multi-criteria 
minimum weight matching problem. The results were established by showing that a multi- 
criteria problem admits an FPTAS if the exact version of the single-criterion problem can 
be solved in pseudo-polynomial time. Let 11 be a single-criterion optimization problem 
with instance set / and objective function w. The exact version of 11 is the following 
decision problem: Given an instance x G / and a number G N, does there exist a 
solution y G sol(a;) with w{x, y) = 
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The exact versions of many single-criterion optimization problems are NP-complete 
since knapsack can be reduced to them easily. But this does not rule out the possibility 
of pseudo-polynomial-time algorithms for them. 

Multi-criteria TSP has been investigated by Ehrgott [11] and Angel et al. [1, 2]. 
Ehrgott [11] considered a generalization of Christofides' algorithm for A-STSP. Instead 
of considering Pareto curves, he measured the quality of a solution y for an instance x 
as a norm of the vector {wi{x, y), . . . , Wk{x, y)). Thus, he encoded the different objective 
functions into a single one, which reduces the problem to a single-criterion problem. The 
approximation ratio achieved is between 3/2 and 2, depending on the norm used to com- 
bine the different criteria. However, by encoding all objective functions into a single one, 
we lose the special properties of multi-criteria optimization problems. 

Angel et al. [1] considered two-criteria STSP(1,2). They presented a 3/2-approxi- 
mation algorithm for this problem by using a local search heuristic. Finally, Angel et al. [2] 
generalized these results to /c-criteria STSP(1, 2) by presenting a 2 — ^^-approximation 
for A: > 3. Although for every fixed k, the approximation ratio is below 2, it converges to 
2 as fc increases. Thus, the ratio tends to the trivial ratio of 2, which can be achieved by 
selecting any Hamiltonian cycle. These two are the only papers about the approximability 
of Pareto curves of multi-criteria TSP we are aware of. 

1.3 Our Results 

All our results hold for an arbitrary but fixed number of objective functions. 

We present a deterministic polynomial-time algorithm that computes (2 -|- e)-approxi- 
mate Pareto curves for A-STSP (Section 2.1). This is the first efficient algorithm for 
computing approximate Pareto curves for this problem. In fact, we show the following 
more general result: If the edge weights satisfy 7-strengthcned triangle inequality for 
7 G [|, 1], then the algorithm computes a (min{l + 7, 2^2^'^^_)_i } + e)-approximate Pareto 
curve for arbitrarily small e > in polynomial time. 

We generalize Christofides' algorithm [10] (cf. Vazirani [26, Sect. 3.2]) to obtain a 
randomized approximation algorithm for multi-criteria A(7)-STSP (Section 2.2). For 
7 G [|)1]) our algorithm achieves an approximation performance of g^l J^^'^^^ -|- e. For 
7 = 1, this yields a ratio of 2 -|- e. 

We consider cycle covers in Section 3. Cycle covers play an important role in the 
design of approximation algorithms for the TSP. We prove that there exists an FPRAS 
for computing approximate Pareto curves of multi-criteria cycle covers. Subsequently, we 
extend this result and show that the multi-criteria variant of the problem of finding graph 
factors of minimum weight admits an FPRAS, too. 

Finally, we analyze a randomized cycle-cover-based algorithm for multi-criteria TSP 
(Section 4): We start by computing an approximate Pareto curve of cycle covers. Then, 
for every cycle cover in the set computed, we remove one edge of every cycle and join 
the paths thus obtained to a Hamiltonian cycle. We analyze the approximation ratio of 
this algorithm for A(7)-STSP (Section 4.2, approximation ratio i^^^2i^^2 + £ for 7 < 1), 

A(7)-ATSP (Section 4.3, ratio \ + + e for 7 < 1/^3), STSP(1,2), and ATSP(1,2) 

(Section 4.4, ratios 4/3 and 3/2, respectively). 

As far as we know, our algorithms are the first approximation algorithms for Pareto 
curves for A-STSP, A(7)-STSP, A(7)-ATSP, and ATSP(1,2). Furthermore, we achieve 
a better approximation ratio for STSP(1,2) than the approximation algorithms by Angel 
et al. [1,2] for aU k. 
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Algorithm 1 The tree doubling algoritiim for multi-criteria A-STSP. 

Input: undirected complete graph G = iV^E); k edge weight functions wi : E ^ N 

(i G [k\)\ £ > 

Output: an approximate Pareto curve T'xsp multi-criteria STSP 

1: compute a (1 + |)-approximate Pareto curve T^^fgT MST on G using the algorithm 

by Papadimitriou and Yannakakis [23] 
2: for all trees T G "Pmst 

3: duplicate all edges in T to obtain an Eulerian graph T 
4: obtain a Hamiltonian cycle S from T by taking shortcuts 
5: put S into 'P:psp 
6: end for 



2 Metric TSP 

In this section, we present two algorithms for A-STSP and A(7)-STSP. Another ap- 
proximation algorithm that can be used for approximating A(7) -STSP, which is based on 
computing cycle covers, will be presented in Section 4. 

The analyses of the algorithms in this section exploit the following result due to 
Bockenhauer et al. [8]. 

Lemma 2.1 (Bockenhauer et al. [8]). Let G = {V, E) be an undirected complete graph with 
an edge weight function w satisfying "/-strengthened triangle inequality for some 7 G 1). 

Let w„ia.x = niaxeg£'(it)(e)) and Wmin = niineg_B(i(;(e)) be the weights of a heaviest and 
lightest edge, respectively. Then ^"^^ < 

Let e and e' be two edges with a common endpoint. Then < . 

Furthermore, we observe the following: Omitting two edges by taking a shortcut re- 
duces the weight by at least 2 • (1 — 7) • Wmin'- The reason is that the two edges {u,v) 
and (w, x) are replaced by (u, x) and w{u, x) < 7 • {w{u, v) + w{v, x)). Thus, the weight is 
leducedhy atleast w{u,v)+w{v,x)—w{u,x) > {l—j)-{w{u,v)+w{v,x)) > 2-{l—j)-Wtam- 



2.1 The Generalized Tree Doubling Algorithm 

Consider the following approximation algorithm for single-criterion A-STSP, which was 
first analyzed by Rosenkrantz et al. [24] (cf. Vazirani [26, Sect. 3.2]): First, we compute a 
minimum spanning tree. Then we duplicate each edge. The result is an Eulerian graph. 
We obtain a Hamiltonian cycle from this graph by walking along an Eulerian cycle. If we 
come back to a vertex that we have already visited, we omit it and take a short-cut to the 
next vertex in the Eulerian cycle. In this way, we obtain an approximation ratio of 2 for 
single-criterion A-STSP. Algorithm 1 is an adaptation of this algorithm to multi-criteria 
STSP. In the following, we estimate the approximation performance of this algorithm. 

Theorem 2.2. For all 7 G [^, 1], Algorithm 1 computes a (min{l -|- 7, 2^2^^27+1 } 
approximate Pareto curve for multi- criteria A (7) -STSP in time polynomial in the input 
size and \/e. 

Proof. We present two analyses showing approximation ratios of l-|-7-|-£ and i2q^2y+i'^^'> 
respectively. The first analysis holds for 7 G [^,1] while the second one only holds for 
7 G [\, 1). However, 1-^7= ^-y-^^-L+i ^r 7 = 1. 
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The key observation of the first analysis is the following: Let T G T^mst' ^^"^ ^ ^® 
any edge in T. Then e appears twice in T, but e cannot appear twice in S since S is a 
Hamiltonian cycle. (We assume that G contains at least three vertices.) Thus, at least 
one copy of e is omitted. This is the moment at which the strengthened triangle inequality 
comes into play. Let ei, 62, . . . , with Cj = {t'j-i, Vj} be a path along the Eulerian cycle 
in T such that this path is replaced by the edge {vq,v^} by taking a shortcut. Then we 
have 

Wi{{vQ, V(}) < 7 ■ {wi{ei) + Wi{e2) + • . . + Wi{et)) 

by iteratively applying 7-strcngthened triangle inequality. (We exploit the fact that 7*^ < 7 
for all c > 1.) Overall, every edge that we omit contributes at most a fraction of 7 of its 
weight. Since we omit at least one copy of every edge e, the two copies of e contribute at 
most (1 + 7) • Wi{e) to S. Thus, 

Wi{S)<{l + -f)-Wi{T) 

for all i e [k]. 

To estimate the overall approximation performance, let S' be an arbitrary Hamiltonian 
cycle. By omitting one edge, we obtain a tree T' . Since P^|sT is a (1 + e/2)-approximate 
Pareto curve for multi-criteria minimum-weight spanning trees on G, there exists a tree 
T G Pm^st with 

w.{T)<{l + '-).w.{T')<{l + '-).w,{S') 
for all i G [k]. Let S be the Hamiltonian cycle obtained from T, then 



'MS) < (1 + 7) • Wi{T) < (1 + 7 + £) . Wi{S') 



for all i e [k]. 

For the second analysis, let again S' be an arbitrary Hamiltonian cycle. This analysis 
only holds for 7 < 1 since it'max/it'min can be unbounded for 7=1. All arguments 
hold simultaneously for all criteria i G [k]. Without loss of generality, we assume that 
mine6Ew;j(e) = 1 for all i E [k], i.e., Wmin = 1- By removing one edge of <S", we obtain 
a tree with a weight of at most Wi{S') — Wmm = Wi{S') — 1. Thus, there exists a tree 
T G 'Pmst from which we obtain a Eulerian graph T with 

Wiif) < 2 • (1 + |) • {w,iS') - 1) = (2 + e) • {w^iS') - 1). 

Let n = \V\ be the number of vertices of the whole graph. Then T contains 2n — 2 edges. 
Thus, in order to obtain a Hamiltonian cycle, we have to remove n — 2 edges by taking 
shortcuts. Every shortcut decreases the weight by at least 2 • (1 — 7). Hence, 

Wi{S) < «;i(f)- (n-2) •2- (I-7) 

< (2 + e) • Wi{S') - (2 + e) - (n - 2) • 2 • (1 - 7) 

< (2 + e)-«;i(5')-2n-(l-7)-(47-2) 

< {2 + e)-Wi{S')-2n-{l-j) 



since 27 > 1. We have Wi{S') = ^l!^y^ for some appropriately chosen a > 1, which 
implies 

w,{S) _ 2n-(l-7)'-a ^ , (1 - 7)' • a 



<2 + £ = 2 + £- 

Wi{S') 27^n 7 



2 
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Algorithm 2 A generalization of Christofides' algorithm for multi-criteria A-STSP. 
Input: undirected complete graph G = iV^E); k edge weight functions wi : E ^ N 

(i G [k\)\ £ > 

Output: an approximate Pareto curve "P^SP *° multi-criteria STSP (with a proba- 
bility of at least 1/2) 

1: compute a (1 + |)-approximate Pareto curve 'P^grp for MST on G using the algorithm 
by Papadimitriou and Yannakakis [23] 
let p be the number of trees in P^g-p 
for all trees T G Pmst 

let V^dd C y be the set of vertices of odd degree in T 

compute PMatch(-^) such that PMatchl-^) ^s a (l + I) -approximate Pareto curve for 
the minimum-weight matching problem on the complete graph induced by V^dd 
with a probability of at least 1 ~ ^ using the algorithm by Papadimitriou and 
Yannakakis [23] 
for all matchings M G V^l^^^{T) do 

let 5 be a Hamiltonian cycle obtained from T U M by taking shortcuts 
put S into T^TSP 
end for 
end for 
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Since Wi(S) < nj^, we also have ^'^g,-^ < a. Thus, the approximation ratio achieved is 



max min q, 2 + e — „ 

a>l \ \ T 

/ / (l-7)2.a\ , 
< max ( mm I a, 2 ^ ) ) + ^ 

= .2 1 . . +^^ 



- 27 -Fl 

which completes the proof of the theorem. □ 

For small values of 7, the bound of 2-^2^^^^ -|- £ is the stronger one, while 1 -|- 7 -|- e 
yields a better bound in case of 7 > I/V2. 

Corollary 2.3. Algorithm 1 computes {2+ e)- approximate Pareto curves for multi-criteria 
A-STSP in time polynomial in the input size and 1/e. 



2.2 A Generalization of Christofides' Algorithm 

In this section, we generalize Christofides' algorithm to multi-criteria A-STSP, which is 
the best approximation algorithm for single-criterion A-STSP known so far. This al- 
gorithm computes approximate Pareto curves of matchings. In case of single-criterion 
A-STSP, we can always find a matching with a weight of at most half of the weight of 
the optimal Hamiltonian cycle. This is in contrast to multi-criteria A-STSP, where the 
weights of the matchings can be arbitrarily close to the weight of the optimal Hamiltonian 
cycle. The reason is that we cannot choose the lighter of two different matchings since 
multiple objective functions are involved; the term '"lighter"' is not well defined. There- 
fore, we only get an approximation ratio of roughly two in this case. But for A(7) -STSP, 
we can show a better upper bound. 
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Theorem 2.4. For 7 G [5, 1], Algorithm 2 is a randomized {-^rz^^ + s) -approximation 
algorithm for multi- criteria A(7)-STSP. Its running time is polynomial in the input size 
and 1/s. 

Proof. The proof consists of two parts. First, we estimate the approximation performance, 
given that all Pareto curves computed arc (H-e/2)-approximate Pareto curves. Second, we 
estimate the success probability, i. e., the probability that such a Pareto curve is computed. 

We assume that all Pareto curves that have to be computed during the execution of 
the algorithm were computed successfully, i.e., with an appropriate approximation ratio. 
Let S' be an arbitrary Hamiltonian cycle of G. Without loss of generality, we assume 
again that Wmm = 1- All the arguments in the following hold for all i. 

There exists a tree T G P^^T ^^^^ 

< (1 + 1) • K(50 - 1). 

Let V^dd be the set of vertices of odd degree in T, and let nodd be its cardinality. Note that 
nodd is even, thus perfect matchings exist on the complete graph induced by V^dd- Let 
n = \ V\ be the number of vertices of the whole graph. Let 5odd be the Hamiltonian cycle 
obtained from S' by taking shortcuts. We get two matchings Mi and M2 on Fodd from 
'S'odd by putting the edges of S'odd alternately into M\ and M2. By 7-triangle inequality, 
we have 

Wi{Mi) + Wi{M2) = u)i(S'odd) < Wi{S') - (n - nodd) • 2 • (1 - 7) 

since every shortcut reduces the weight by at least (1 — 7) • 2i(;min- Now, Wi{Mi) • < 
Wi{M2) according to Lemma 2.1. (Note that we do not make any assumptions whether 
Ml or M2 is the lighter matching.) There exists a matching M G ^Match(^) ^^^^ 

Wi{M)<{l + ^)-Wi{Mi). 

Hence, 

^^ = Wi{M).{l + ^) < {l+'-).{w,{M,) + w,{M2)) 

< (1 + |) • {wi{S') - 2 • (n - nodd) • (1 - 7)) • 

Let D be Eulerian graph obtained by taking the union of the tree T and the matching M. 
By the arguments above, we can bound its weight as follows: 

Wi{D) = Wi{T)+Wi{M) 

< (1 + |) • {wi{S') - 1 + 7 • {wiiS') - 2 • (n - nodd) • (1 - 7))) 
= (1 + |) • ((1 + 7) • w,iS') - (1 + 27 • (n - nodd) • (1 - 7)))- 

The Eulerian graph D consists of n — 1 + nodd/2 edges, the tour S constructed from D 
consists only of n edges. Thus, nodd/2 — 1 edges are removed and, by 7-triangle inequality, 
we have 

Wi{S) < Wi{D) - 2 • (nodd/2 - 1) • (1 - 7) = Wi{D) - (1 - 7) • ("odd - 2). 
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By combining these inequalities, we obtain 



w. 



{S) < {l+'-)-{{l + j)-w,{S') 

-(1 + 27 • (n - nodd) • (1 - 7) + (1 - 7) • Kdd - 2))) 
= (l + |)-((l + 7)-^.(5') 

( 1-2^+ 27 +(l - 7)-(^-(n - nodd) + "odd))) 

>o >i 

< (l + £).((l+7).^,(5')-n.(l-7)). 



Now we have t^max < since i«min = 1- Thus, Wi{S') < ■ n. We choose a > 1 such 

(l-7)-a' 



that Wi{S') = jTz\~^ which imphes 



1-7 

The last inequality holds since (1 + 7) - ^%=?^ < 1 + 7 < 2. 



27^ 

Since Wi{S) < nj^, we also have < a. Thus, the approximation ratio achieved 



IS 



«(l- 7)' 

a>l" V V""" ' " 272 



max ( min ( a, (1 + 7) h e 



^ 27^ + 27^ 27^ ^ 27^ + 27^ 

- 372-27 + 1 372-27 + 1 - 372-27 + 1 

We obtain the first inequality by observing that (1 + 7) — "^^"^^^ — h £ is monotonically 
decreasing in a: The maximum of the minimum is therefore assumed for a = (1 + 7) — 

°'^^2^2^ — h £. The second inequality follows from the fact that 3^:^^^^_|_x ^ 1 for 76 [^5 !]• 
The analysis so far holds only for 7 < 1 since for 7 = 1, division by zero occurs at 
some points in the analysis. For 7 = 1, we obtain a ratio of 2 + e = ^^'^ J^^"^-^ + e: We have 
w{M) < (1 + f ) • w{S') and w{T) < (l + |) • w{S'), which implies the bound. 

What remains to be proved is that the algorithm succeeds with a probability of at least 
1/2. First, we observe that if we iterate the randomized computation of an approximate 
Pareto curve, then we do not have to decide which set is indeed such an approximate 
Pareto curve. Instead, we can take the union of all solutions computed and remove all 
dominated solutions of the set thus obtained. The only randomization in Algorithm 2 is 
in the computation of the approximate Pareto curves for the matching problems. The 
number of curves to be computed is p, which is bounded by a polynomial of the input 
size and 1/e. We can achieve a failure probability of at most ^ by performing log(2p) 
iterations of the FPRAS for the matching problem. The probability that one of the Pareto 
curve computations fails is thus at most p ■ ^ = 1/2, which completes the proof of the 
theorem. □ 

We compare the ratios obtained by the two algorithms of this sections and the cycle 
cover algorithm of Section 4 in Section 5.1. 
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3 Matchings and Cycle Covers 



A cycle cover of a graph is a spanning subgraph that consists solely of cycles such that 
every vertex is part of exactly one cycle. Many approximation algorithms for the single- 
criterion TSP are based on cycle covers. These approximation algorithms usually start by 
computing an initial cycle cover and then join the cycles to obtain a Hamiltonian cycle. 
This technique is called subtour patching [15]. We show that there exist FPRASs for 
multi-criteria cycle cover problems. 

ACC, the cycle cover problem in directed graphs, is equivalent to finding matchings 
of minimum weight in bipartite graphs (assignment problem). An FPRAS for the multi- 
criteria matching problem is also an FPRAS for the multi-criteria matching problem in 
bipartite graphs. Hence, multi-criteria ACC also admits an FPRAS. 

Theorem 3.1. There exists an FPRAS for multi- criteria ACC. 

To show that multi-criteria SCC admits an FPRAS, we show that arbitrary graph 
factor problems admit FPRASs. Let G = (V, E) be a graph and / : F — >^ N be function. 
A subset F C is called an /-factor of G if all vertices v &V have a degree of exactly 
f{v) in the graph iy,F). 

Cycle covers of undirected graphs are also known as two-factors since every vertex is 
incident to exactly two edges. Thus, they are a special case of graph factors. 

The graph factor problem GFP is the following minimization problem: An instance is 
an undirected graph G = {V, E) with a function / : F ^ N and an edge weight function 
w : E ^ 'H. The aim is to find an /-factor of minimum weight. 

To show that multi-criteria GFP, and thus multi-criteria SCC as well, admits an 
FPRAS, we exploit Tutte's reduction [25], which reduces arbitrary graph factor problems 
to matchings (matchings arc also known as one-factors since every vertex is incident to 
exactly one edge of the matching). We omit a description of the reduction, but refer to 
Lovasz and Plummer [19] or Tutte [25] for the details. Overall, we obtain the following 
result. 

Theorem 3.2. Multi- criteria GFP and multi- criteria SCC admit an FPRAS. 



4 Approximations Based on Cycle Covers 
4.1 The Algorithm 

The generic outline of a cycle-cover-based algorithm is the following: Start by computing 
a cycle cover. Then remove one edge of every cycle. Finally, join the paths thus obtained 
to form a Hamiltonian cycle. 

Algorithm 3 is our generalization of this algorithm to multi-criteria TSP. It achieves 
a constant approximation ratio if the quotient of the weight of the heaviest edge and the 
weight of the lightest edge is bounded. 

In this section, we present a general analysis of the approximation ratio of this algo- 
rithm. We will refine the analysis for multi-criteria A(7)-STSP (Section 4.2) to get an 
improved approximation ratio. Furthermore, we apply the analysis to get approximation 
results for multi-criteria A(7) -ATSP (Section 4.3) and STSP(1, 2) and ATSP(1, 2) (Sec- 
tion 4.4). We analyze Algorithm 3 in terms of the number an of edges that have to be 
removed and the quotient /? = Wmax/'U^min- 

Lemma 4.1. Assume that at most an edges have to be removed from each cycle cover 
and that ZZ^'^'^Tm < P for all ie[k]. 
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Algorithm 3 An approximation algorithm for multi-criteria TSP based on cycle covers. 
Input: complete graph G = {V^E)] k edge weight functions Wi {i € [k]); e' > 
Output: an approximate Pareto curve P^sp to multi-criteria TSP (with a probability of 
at least 1/2) 

1: compute a (1 -|- e')-approximate Pareto curve Vqc to the multi-criteria cycle cover 

problem on G using the algorithm by Papadimitriou and Yannakakis [23] 
2: for all cycle covers C G Vcc do 
3: for all cycles c of C do 
4: remove one edge of c 

5: end for 

6: join the paths to form a Hamiltonian cycle S 
7: add S to V^fp 
8: end for 



Then Algorithm 3 is a randomized (l -|-a(/3 — 1) -|-e) approximation algorithm for every 
e > 0. Its running-time is polynomial in the input size and 1/e. 

Proof. Without loss of generality, let miueg e Wi{e) = 1 for all i G [k]. We run the algorithm 
with some e' that depends on a, /3, and e and that we will specify later on. Let S' be 
an arbitrary Hamiltonian cycle. Then there exists a cycle cover C in Vcc with w{C) < 
(1 -I- e') ■ w{S'). We obtain a Hamiltonian cycle S from C such that 

WiiS) < WiiG) + aniP - 1) 

for all i £ [k]. The reason for this is that every edge removed has a weight of at least 1 
and every edge added has a weight of at most /3. Now we have for all i G [k] 



< (l + eO- " + ""f "'^ = (l + .0-(l + a(/3-l)) 

< l + a{(3-l)+s 

for e' < jqi^^pujy, which proves the lemma. □ 



4.2 Refined Analysis for A(7) -STSP 

Prom the general analysis (Lemma 4.1), we obtain an approximation ratio of | -|- | • 

2 

-f^ + e for A(7) -STSP. In this section, we present a refined analysis that yields a better 
approximation ratio. 

Consider any cycle c of a cycle cover of Vcc- Then there will be an edge sr of c 
that will be removed and an edge adjacent to that will be added during the joining 
process. Finally, there exists an edge of c that is adjacent to both and (Figure 1 
shows an example). Note that while cr is uniquely determined, once the edges have been 
removed and the new edges have been added, the edge is not since there are two edges 
that connect c to other cycles of the cycle cover. However, once we have fixed ca for one 
cycle c, the corresponding ex is uniquely determined, and the e'^ and e'^ of all other cycles 
c' are also determined. 

By Lemma 2.1, we have Wi{eR) > ^-^-Wi^eA) and Wi^ex) > ^-^■Wi{eA)- AH arguments 
in the following hold for all weight functions simultaneously. Thus, we restrict ourselves 
to considering one fixed weight function Wi for some i G [k] to simplify the arguments. 
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c 



eA 



eR 



c 



(a) Cycle cover, before the patching. 



(b) Hamiltonian cycle, after the patching. 



Figure 1: Two cycles c and c' before and after joining the cycles to a Hamiltonian cycle. 
The edges eji, ex, and belong to c while e'^, e'^, and e'^ belong to d . 



Let wr be the total weight of edges removed, wa be the total weight of edges added, 
and wk be the total weight of edges of C and S that are adjacent to edges added. Then 
we have w{C) = w + wk + wr for some suitably chosen w >Q, which is the total weight 
of all edges not taken into account so fax. Thus, 

Wi{S) W + WK + WA , , WA-WR 

TTTT ~ = J- H = tt. 

Wi(L) w + Wk + Wr w + wk + wr 

Since R is monotonically decreasing with respect to wr, we obtain 

WA-^ - WA 
R<1 + \- = B!. 

W + WK + ^ ■ WA 

Exploiting further that E! is monotonically decreasing in wk-, we get 

w + ^ ■ 2wA jw + [l- j)2wA 

The inequalities wa < 3^1-^) ™d w > n/3 hold since every cycle has a length of at least 
three. We exploit the fact that R" is monotonically increasing with respect to wa and 
monotonically decreasing with respect to w: 

^.^,1 3g^-(^7-l) ^, |g^-(27-l) 

^ + (1-7) •2- 3^ 7 + 47^ 1 + 37-472 

We run the algorithm with some e' > that depends on 7. We will specify e' in a 
moment. Let S' be an arbitrary Hamiltonian cycle and C G Vcc be a cycle cover with 
Wi{C) < (1 + e') ■ w{S') for all i G [k]. Let S be the Hamiltonian cycle obtained from C. 
Then 

w,{S)<{l + e'). ^^\^\^^ .w.{S'). 

For a given e > 0, we choose e' such that jqij^rj^ < £• The set "P^sp is a ( jqij^rj^+s)- 
approximate Pareto curve with a probability of at least 1/2, which implies the following 
theorem. 

Theorem 4.2. For 7 G [\, 1), Algorithm 3 is a randomized {j^^z^ -approximation 
algorithm for all £ > 0. Its running-time is polynomial in the input size and 1/e. 

In Section 5.1, wc compare the approximation ratios of the cycle cover algorithm for 
A(7)-STSP to the tree doubling and Christofides' algorithm. 
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Figure 2: The approximation ratio of Algorithm 3 achieved for A(7) -ATSP subject to 7 
compared to the trivial approximation ratio of i«max/i«min- 

4.3 The Cycle Cover Algorithm for A(7)-ATSP 

For multi-criteria A(7) -ATSP, our algorithm yields a constant factor approximation if 
7 < since it'max/^^min is bounded from above by -^z^^ for such 7. For larger values of 
7, this ratio can be unbounded. 

Lemma 4.3 (Chandran and Ram [9]). Let 7 G [1/2, 1). Let G = {V, E) he a directed com- 
plete graph, and let w : E ^ 'H he an edge weight function satisfying '^-triangle inequality. 
Let u;min = mineeE w{e) and lUmax = maxggB w{e). 

If^< l/-\/3, then < t^-^. If J > l/-\/3, then can be unhounded. 

By combining Lemma 4.1 and Lemma 4.3, we obtain the following result. 

Theorem 4.4. For 7 < l/\/3; Algorithm 3 is a randomized (^+ ^^^3^2 +s)- approximation 
algorithm for A(7)-ATSP. Its running-time is polynomial in the input size and 1/e. 

Figure 2 shows the approximation ratio achieved for multi-criteria A (7) -ATSP subject 
to 7 and compared to the trivial ratio of u'max/ii'max- 

We leave as an open problem to generalize the analysis to larger values of 7. However, 
it seems to be hard to find a constant factor approximation for 7 = 1, i.e., for multi- 
criteria A -ATSP, since this would immediately yield a constant factor approximation for 
single-criterion A -ATSP. 

4.4 TSP with Weights One and Two 

Now we analyze the cycle cover algorithm for multi-criteria TSP with weights one and two. 
For both STSP(1,2) and ATSP(1,2), we have (3 = 2, i.e., Wmax/^min = 2. Furthermore, 
for STSP(1, 2), we have a < 1/3, while we only have a < 1/2 in case of ATSP(1, 2). The 
approximation ratio follows by exploiting Lemma 4.1. 

Note that the edge weights and thus the objective functions are polynomially bounded 
for STSP(1,2) and ATSP(1,2). Thus, wc can compute a Parcto curve of cycle covers 
instead of only a (1 + e)-approximatc Parcto curve. This implies that we do not have an 
additional e in the approximation ratios in the following theorems. 

Theorem 4.5. Algorithm 3 is a randomized A /'i- approximation algorithm for multi-crite- 
ria STSP(1,2). Its running-time is polynomial. 




cycle cover 
trivial ratio 
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Figure 3: Approximation ratios subject to 7 achieved by the tree doubhng algorithm 
(Algorithm 1), Christofidcs' algorithm (Algorithm 2), and the cycle cover algorithm (Al- 
gorithm 3, Section 4), for which both the ratio obtained from the general analysis (Sec- 
tion 4.1) and from the refined analysis (Section 4.2) are shown. 

Theorem 4.6. Algorithm 3 is a randomized S/2-approximation algorithm for multi-crite- 
ria ATSP(1,2). Its running-time is polynomial. 

5 Concluding Remarks 

5.1 Compciring the Approximation Ratios 

Let us compare the approximation ratios for A(7)-STSP achieved by the tree doubling 
algorithm (Algorithm 1), Christofides' algorithm (Algorithm 2), and the cycle cover algo- 
rithm (Algorithm 3). 

Figure 3 shows the approximation ratios achieved by these algorithms subject to 7. 
Figure 4 shows the approximation ratios achieved deterministically (by the tree doubling 
algorithm) and randomized (by a combination of Christofides' and the cycle cover algo- 
rithm). The ratios are compared to the trivial ratio of ^t;max/^^min and to the currently 
best known approximation ratio for single-criterion A(7)-STSP. Note that in particular 
for small values of 7, our algorithms for multi-criteria A(7) -STSP come close to achieving 
the ratio of the best algorithms for single-criterion A (7) -STSP. 

5.2 Open Problems 

Our approximation algorithm for multi-criteria A(7)-ATSP works only for 7 < l/\/3- 
Thus, wc arc interested in finding constant factor approximation algorithms also for 7 > 
l/\/3, which exist for all 7 < 1 for single-criterion A(7)-ATSP [7,9] 

The cycle-cover-based algorithm for Max-TSP, where Hamiltonian cycles of max;imum 
weight are sought, does not seem to perform well for multi-criteria Max-TSP. The reason 
for this is that the approximation algorithms for Max-TSP that base on cycle covers usually 
contain a statement like "remove the lightest edge of every cycle" . While this works for 
single-criterion TSP, the term "lightest edge" is not well-defined for multi-criteria traveling 
salesman problems. We are particularly curious about the approximability of multi-criteria 
Max-TSP. 
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Figure 4: Approximation ratios subject to 7. The deterministic ratio is achieved by the 
tree doubling algorithm. Combining Christofides' and the cycle cover algorithm yields the 
randomized ratio. For comparison, the current ratio for single-criterion A(7) -STSP and 
the trivial ratio Wmax/wmin are also shown. 
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